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We extend to the spectral function an approach which allowed us to calculate the quasiparticle 
weight for destruction of a real electron Z ccr (k) (in contrast to that of creation of a spinless holon 
Zh(k)) in a generalized t — J model, using the self-consistent Born approximation (SCBA). We 
compare our results with those obtained using the alternative approach of Sushkov et al, which 
also uses the SCBA. The results for Z C(T (k) are also compared with results obtained using the string 
picture and with exact diagonalizations of a 32-site square cluster. While on a qualitative level, all 
results look similar, our SCBA approach seems to compare better with the ED one. The effect of 
hopping beyond nearest neighbors, and that of the three-site term are discussed. 



I. INTRODUCTION 

The momentum distribution function of holes in a 
quantum antiferromagnetic background has been a sub- 
ject of considerable interest since the discovery of high-T c 
systems. Reliable information on these quantities was ob- 
tained from the fix#ct diagonalization of the t — J model 
in small clustersEJu. An analytical approximation which 
brings considerable insight into the underlying physics is 
based on the string pictured. Within the string picture, 
for realistic t > J, the movement of the hole can be sep- 
arated into a fast motion around a fixed position j on 
the lattice (to which the hole is attracted by a string po- 
tential caused by the distortion of the Neel background), 
and a slow motion of j due essentially to spin fluctuations 
which restore the Neel background as j is displaced. The 
resulting quasiparticle weight as a function of wave vector 
Z ca (k) for one hole, agrees very_well with exact diagonal- 
izations (ED) of a 4 x 4 clustem. However, this cluster is 
still too small and finite-size effects are importantoEl. 

Another successful analytical approaQh-JS the self- 
consistent Born approximation (SCBA)tjlHl. The re- 
sulting dispersion of one hole in the t — J model is 
in very good|-aepeement with exact diagonalizations of 
small clustersBEHlTJ, including a square cluster of 32 sites, 
which contains the most important symmetry points and 
smaller finite-size effects in comparison with previous 
calculationsQ. The SCBA has also been used to cal- 
culate the doping dependence of the superconducting 
critical temperature T c in qualitative agreement with 
experiment!^. However, due to the particular representa- 
tion used (see Eq. (|J) below), the Green function which 
results from the SCBA is that of a spinless holon Gh, 
while the real particle becomes related with a composite 
operator Cka (composed of a holon and a spin deviation) . 
The holon quasiparticle weights Zh(k) differ from Z ca {k) 
which are the physical quantities calculated by ED and 
accessible to experiment. Only recently, motivated hy 
photoemission experiments in insulating Sr 2 Cu02Cl2tJ, 



two approaches appeared which relate, Z ca and the 
Green function of the physical hole G C(T E3 with Gh within 
the framework of the SCBA. In addition, to fit accurately 
the experimentally observed dispersion, it is necessary to 
include second- (£2) and ..third- nearest-neighbor hopping 
(t 3 ) to the t — J modeOO, and to explain qualitatively 
the observed intensities, it is necessary to-.pansider the 
strong-coupling limit of a Hubbard modelli-ffiS. This im- 
plies that a three-site term t" = — J/4 should be included 
in the model, and the one-particle operators should be 
transformed. While a generalized Hubbard model is able 
to provide a consistent picture of the-pbserved charge and 
and spin excitations in Sr 2 Cu0 2 Cl 2 EJ3, we must warn that 
the effective strong-coupling low-energy .effective model 
derived from a realistic multiband mode£3~E3, although 
also includes i 2 and £3, can have a different t", even of 
opposite sign, favoring instead of suppressing of d-wave 
superconductivity and a resonance-valence-bond ground 
stateMtll. 

For future theoretical .studies, as well as to compare 
with photoemissionll3c3'Ej, or Qthp: experiments (like 
electron-energy-loss spectroscopyE3EII) it is important to 
compare the four above mentioned approaches to calcu- 
late the quasiparticle dispersion and intensities, trying to 
establish their relative accuracy or convenience. This is 
the main purpose of this work. 

In Section 2 we derive an expression for the spectral 
density of the physical hole p ca = —\a\G ca l^ in terms 
of the Green function of the spinless hole Gh within the 
SCBA, and briefly describe the alternative expression de- 
rived by Sushkov et a/.lla for G ca . Section 3 contains a 
comparison of both resulting p C(J (k,Ci;) after solving the 
SCBA equations. In Section 4, we compare the results 
for Z ccr (k) obtained using both SCBA approaches with 
ED results for the t — J modelO, and those derived from 
a string pictureQ. In Section 5 we include £2, £3 and £1 
and compare recent ED results for the 32-site cluster^ 
with the SCBA ones. Section 6 contains the conclusions. 
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II. GREEN FUNCTION OF A PHYSICAL HOLE 
WITHIN THE SCBA 

We consider the dynamics of one hole in a square lat- 
tice described by a generalized t — J model: 



H =-J2 ^4+8*°™ + g E 
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The first term contains hopping to first, second and third 
nearest neighbors with parameters t\, t%, £3 respectively. 
The nearest neighbors of site i are labeled as i + r\. In the 
SCBA, long-range antiferromagnetic order is assumed, 
and a spinless fermion hj at each site is introducedLTHl. 
Calling A (B) the sublattice of positive (negative) spin 
projections in the Neel state, one can use the following 
representation: 



Cjf = hj, ai = h\a t , if i e A 
Cj T = h\ai, Cii = h\ , if i e B, 



(2) 



where a\ creates a spin deviation at site i in the Neel 
state, and there is a constraint that at the same site 
there cannot be both, 2. hole and a spin deviation. This 
constraint is neglected!! The exchange part of Eq. (1) 
is diagonalized by a standard canonical transformation, 
and retaining only linear terms in the spin deviations for 
the other terms, the Hamiltonian takes the form: 



H =} j u q a\a q + ^ e kh{h k 
9 



-|= Af (k, q) {h\h k - q a q + H.c), 



(3) 



kq 



where a q = u q a q - v q a q , uj q = 2 J v q , e k = (t 2 + 2(1 - 
x)t")e 2 {k) + (t 3 + (l-x)t")e 1 (2k) x= 1/iVis the doping, 
ei(fc) = 47 fc , e 2 (k) = 4 cos k x cos k y , 



v 2 q = a 



2 - 1 = — - - v = (1 -<r 2 ) 1/2 



7 9 = (cosfe + cos%)/2, sgn(w g ) = sgn( 7(3 ), 
M(k, q) = u qlk - q + v q ~/ k , (4) 

and u q > 0. 

From the Hamiltonian Eq.(|J), the SCBA allows to cal- 
culate the holon Green function Gh accurately through 
the self-consistent solution of the following two equations: 



G h 1 (k, w) = a; - efc - S(k, w) + ie. 



= 4 is the coordination number. 



(5) 



However, the physical operator is the Fourier trans- 
form of and since it is a composite operator in 
the representation Eq.(||), it is not trivial to find its 
Green function G ca . Recently, two different .approaches 
to calculate the quasiparticle weight Z C(T (k)li-J and the 
whole Green function G ccr (k, uj)i3 using the SCBA have 
been proposed. Here we extend to the spectral density 
p C(T (k,w) = — ImG C(T (k,Ci;)/7r our previous derivatiorO. 
We assume for the moment that the system is finite and 
its eigenvalues are discrete. The idea is to relate G ca and 
Gh , using the equation of motiorH to find the wave func- 
tion for each eigenvector. If is an eigenstate of Eq. 
(^) with total wave vector k and other quantum numbers 
labeled by n, it can be expanded as: 

|V4 l )=^(k)^|0) + ^V^(k,q)/ l t 4| ) 



-l^^(k, qi ,q 2 )^ 
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and the Schrodinger equation (AJJ — H)\ip2) — 0, leads to 
an infinite set of equations for the, A™. Neglecting some 
terms not described by the SCBAEj, the first two of these 
equations are: 



A£(k)(A£ 
A?(k,q)(AJ- 



e k 



- £k-q 

-y 
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-o; 9 )=ijrti[4f(k)M(k,q) 
A^(k,q,q 2 )M(k-q,q 2 )]. (7) 



A solution of the set of approximate equations can be 
solved relating A™ with A^Ea. In particular, if 



^(k, q, q 2 ) = *M?(k, q)M(k - q, q 2 ) 
xG h (k-q-q 2 ,A^-w, 



(8) 



with G/j(k, uj) and S(k, to) satisfying Eqs. (||), from Eqs. 
(0) it follows that: 

A"(k,q) = zt 1 ^ l (k)Af(k,q)G h (k-q,A[!-^), (9) 

and the eigenvalue equation: 

A£(k)[A£-e fe -E(k,A£)] = 0. (10) 

Since we are considering the case of only one hole 
at zero temperature (otherwise Eqs. (||) have to be 
generalizecoEj) , the relevant Fock space is composed by 
|0) (the ground state of the Heisenberg antiferromagnet), 
and the eigenstates with one added hole \ip k ). In this-jje^ 
stricted Hilbert space, the Lehmann representationsOLJ 
of the relevant Green functions read: 



G CCT (k,w) = 



E 



m\ck.m 2 + m +{ ^^\c k .m 



AJ + ie 



K + ie 



(ii) 
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Using Eqs. (0) and the expression of the a% in terms 
of the magnon operators a q to Fourier transform Ci CT E_l, 
Eqs.(||), (||), (|Tl]), and some algebra, we find for the ratio 
of spectral functions: 



(12) 



where the sum runs over all wave vectors of the antifer- 
romagnetic Brillouin zone, excluding q = 0, and v q and 
M (k, q) are given by Eq. (||) . 

Although Eq. ( |l2| ) has been derived assuming that 
the eigenvalues are discrete, we expect it to be also valid 
for a continuum distribution of energy levels (the inco- 
herent background of the spectrum), or in other words, 
when the small imaginary part it in Eqs. (|5|) and ( JT^ ) is 
larger than the average spacing between the levels. This 
seems to be confirmed by the results shown in the next 
section. A similar method of using the Schrodinger equa- 
tion assuming discrete eigenvalues was used to obtain the 
non-interacting Green function of the Anderson modclE3. 

Sushkov et alHH used a different method to relate 
G ca (k, uj) with G/ l (k, uj). They treat the operators Ck a as 
originated by an external perturbation to the Hamilto- 
nian Eq. (^). This perturbation is characterized by two 
vertices in which the physical hole is transformed into a 
spinless hole h\., with or without exchange of a magnon. 
The physical Green function G ca (k,uj) is obtained from 
a Dyson equation. Another difference with our approach 
lies in the different normalization of the operators hi, 
which is equivalent to another way of treating the con- 
straint that there cannot be a hole and a spin deviation 
at the same site. We have neglected it since it has been 
shown that it does not affect the results for a quantum 
antifcrromagnet (the situation is different if the spins are 
described by an Ising model)Q. Sushkov et al. also have a 
factor \[2 in their normalization of the Ck<y , which intro- 
duces a factor 2 in G ca . Dropping this factor to compare 
with other results, their final expression for the physical 
Green function reads: 



G CCT (k, w) = 0.4G h (k, uj) + — 22 v 2 q Gh{^ - q, w - «,) 
+4VoIhG h {k, u) — v q M(k, q)G /l (k - q, u - oj q ) 



III. COMPARISON OF SPECTRAL DENSITIES 
WITHIN SCBA. 

The expression of the spectral density of the physical 
hole p c<T (k-o^j derived from Eq. ( [l3| ) is quite different 
from Eq. (|12]). Nevertheless, in both cases it is necessary 
first to solve numerically the SCBA Eqs. (||) to obtain 
the Green function of the spinless hole G/,. (k, uj) . To per- 
form this task with high accuracy, we have discretized the 
frequencies in intervals of Atu = 10~ 4 <i and have taken 
the small imaginary part as e = 5 Aw. We have chosen a 
cluster of 16 x 16 sites. 
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FIG. 1. Spectral function of the physical hole obtained 
from the SCBA for wave vector k= (37r/4, 7r/4). Full line: 
result of Eq.(12). Dashed line: result derived from the Dyson 
Eq. (13) [13]. Parameters are t\ = 1, J = 0.3, = —0.3, 
tz = 0.2, t" = 0. 




FIG. 2. Same as Fig. 1 for k= (7r/8,7r/8). 
1 ' 

+32tlG h (k,uj)[—^2v q M(k,q)G h (k-(i,uj-uj q )} 2 , (13) In Figs. 1 and 2, we represent the resulting spectral 



where v q and M(k, q) are given by Eq. 



densities for a wave vector k at the boundary of the an- 
tiferromagnetic Brillouin zone, and another one near the 
zone center respectively. Except for wave vectors near 
(■7T, 7r) at frequencies near the quasiparticle pole, the con- 
tribution of the term in t\ in Eq. (12) is small com- 
pared to 1, and as a rough first approximation Eq. ( p^ ) 
gives p ca (k,uj) ~ pi l (k,u;)/2. Thus, the result of our Eq. 



3 



( |l2| ) looks qualitatively-, similar to the known results for 
Ph within the SCBAlEj. There are several peaks which 
can be qualitatively understood within the string pic- 
ture as originated by different bound states of the string 
potential, which acquire dispersion as the center of the 
string potential is displaced by spin fluctuations or terms 
of sixth order in tJB'U. The lowest peak corresponds to 
the coherent quasiparticle state. Note that the structure 
in the incoherent background displaying several different 
peaks cannot be resolved if the imaginary part in Eqs.j-dq) 
is as large as that used by Sushkov et al. (e = O.Kjij). 
The same happens for certain wave vectors and parame- 
ters with the quasiparticle peaks, when the quasiparticle 
energy Afe lies too near the incoherent background. 

The result for p ca -(k,u>) obtained using the Dyson Eq. 
( |T3| ) derived by Sushkov et alS3 looks in general, and for 
any reasonable parameters of the generalized t — J model, 
similar to ours, except for two differences: i) the intensity 
is a little bit smaller for the peaks already present in the 
spinlcss holon result ph(k,u}) (see Fif. 1), ii) new peaks 
appear, which apparently do not have a physical meaning 
and seem to be an artifact of the approximations involved 
in the Dyson equation. In particular, except for high 
symmetry points of the Brillouin zone (like (0,0), (0, 7r) 
or (tt,tt)), for wave vectors k such that the quasiparti- 
cle energy Afe is near A( ,o) (the bottom of the electron 
band), a new quasiparticle peak appears. Comparison 
with the position and intensity of the quasiparticle peak 
at k = (37r/4,7r/4) obtained, from exact diagonalization 
of a square 32-site clusterOo, described in the following 
two sections, confirm that this peak is spurious and is 
disregarded in the following. 



IV. QUASIPARTICLE WEIGHTS IN THE T-J 
MODEL 



To obtain the quasiparticle intensity for each wave vec- 
tor of the physical hole Z ca {k) within the SCBA, we have 
fitted the part of the spectral density (described by Eq. 
@ or Eq. Q and Eq. (g)) in the neighborhood of the 
quasiparticle peak by a sum of several Lorentzian func- 



tions. Since the imaginary part e = 5 x 10 ii we have 
taken in the numerical solution of the SCBA Eqs. (|^) 
is very small, we can isolate the quasiparticle peak (its 
width is practically identical to 2e), even in some cases 
where the distance of this peak to the incoherent back- 
ground is smaller than e. Integrating the corresponding 
Lorentzian we obtain Z C(T (k). We have verified that us- 
ing this method, finite-size effects are practically absent 
in our 16 x 16 cluster. The peaks introduced by Eq. 
( |l3| ) which are absent in the spinless hole spectral den- 
sity ph(k,u>) were neglected. 




(0,0) (jt/2,«2) (n,jt) (x,nl2) (xfi) (n/2,0) (0,0)(0,ji) (it/2,ji/2) (n,0) 



FIG. 3. Top: quasiparticle energies obtained with the 
SCBA (solid squares). Bottom: corresponding quasiparti- 
cle weights obtained using Eq.(12) [12] (solid triangles) and 
from Dyson Eq. (13) [13] (solid circles). Open circles de- 
note the corresponding results of an exact diagonalization 
(ED) of a square cluster of 32 sites [3]. Parameters are 
h = l, J = 0.3, t 2 =t 3 = t" = 0. 

In Fig. 3 we compare quasiparticle energies Afe and 
weights Z C(J (k) with those obtained from an exact diag- 
onalization (ED) of a square cluster of 32 sitestl. The 
SCBA results for Afe were already published!!! We in- 
verted the sign of Afe (the electron instead of the hole 
representation is used) in the following, in order to facil- 
itate comparison with previous calculations and experi- 
ment. We also shifted the Afe in order that they coincide 
for the ground-state wave vector k = (7r/2,7r/2), as in 
Rcf.El. The agreement between ED and SCBA results for 
Afe is excellent along the boundary of the antiferromag- 
netic Brillouin zone. However, there are some discrep- 
ancies near the zone center, particularly for k = (0, 0) 
and k = (7r/4, 7t/4). We ascribe this to finite-size ef- 
fects, since in the thermodynamic limit Afe = \k+(Tr,n)> 
due to the folding of the Brillouin zone caused by the 
antifcrromagnetic symmetry breaking. Except for the 
above mentioned two wave vectors, the disagreement be- 
tween the-results for Z C(T (k) of the ED and our SCBA 
approachH (Eq. (p^)) is less than 7%. Using instead 
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the SCBA expression © of Sushkov et aim, we obtain 
a quasiparticle weight which is ~ 20% below our results. 
Previous comparison of ED results for Z ca (k) on a square 
cluster of 20 sitesa and SCBA results on a 20 x 20 cluster- 
using Eq. (|T^), also agreed very well except at k = (0, 0). 

0.20 




FIG. 4. Quasiparticle weights obtained with the SCBA 
using Eq.(12) [12] (solid triangles) and Eq. (13) [13] 
(solid circles). The full thick line corresponds to results of 
the string picture scanned from Ref. 4. Parameters are 
ti = 1, J = 0.4, t 2 = tz = t" = 0. 

Comparison ^tween results of ED for clusters between 
16 and 32 sitesQU, suggest that finite-size effects are very 
important for the 4x4 cluster. In Fig. 4 we com- 
pare the results for Z ccr (k) obtained using the string pic- 
ture, scanned from Ref.B with both SCBA approaches. 
ED results on large enough clusters are not available for 
these parameters. While the three curves look qualita- 
tively similar, it seems that the string picture results un- 
derestimate Z C(T (k), particularly at k = (0,0) and near 
k = (tt, 7t). 

V. QUASIPARTICLE ENERGIES AND WEIGHTS 
IN THE GENERALIZED T-J MODEL 

In this Section, we compare results for A& and Z C(7 (k) 
which we obtained using the SCBA, with correspond- 
ing recent ED results, in which hoppings beyond near- 
est neighbors t% and the three-site term t" were 
included^. A motivation to include these, pterms is 
that the inclusion of t% and £3 is necessarylla-EJ to ex- 
plain the experimentally observed dispersion in insulat- 
ing S^CuC^Cylj. However, to explain qualitatively 
the observed quasiparticle intensities, it is necessary to 
include ai-iaast Hubbard corrections to the relevant 
operatorsEjni3ii3. Their effect can be included in any 
theoretical approach and its main effect is. increase 
the intensities near the Brillouin zone centeroEjiij. 
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FIG. 5. Same as Fig. 3 for the parameters of Fig. 1. The 
ED results were taken from Ref. 25. 

Fig. 5 contains the results for t" — 0. Comparing 
with Fig. 3 one can see that the effect of ti and £3 
is mainly to shift A( w 0) towards the incoherent back- 
ground, and as a consequence, the weight Z ca {-K, 0) is 
strongly reduced. The same happens for neighboring k. 
Also, the agreement between ED and SCBA results for 
Afc near k = (n,0) is not so good as for t 2 = t 3 = 0. 
In spite of this, the agreement between ED and SCBA 
results for Z ccr (k) is still very good except at the points 
k = (0,0), k = (7r/4,7r/4) and k = (3tt/4, 3tt/4), for 
which finite-size effects are present in the ED calcula- 
tions as evidenced by the fact that A& 7^ A fc+ ( 7r >7r ) . In 
Fig. 6, we include the three-site term t" = — J/4, with 
magnitude corresponding to the strong-coupling limit of 
the Hubbard model. Its effect is to lower the quasipar- 
ticle energies Afe near the wave vectors (0,0) and (tt, 0), 
increasing the total dispersion. While the comparison 
between SCBA and ED results for Afc is not affected ap- 
preciably by t" , the corresponding weights Z ca (k) clearly 
disagree for wave vectors (7r/2, 0) and (tt, 7t/2). However, 
the agreement between weights on the boundaries of the 
antiferromagnetic Brillouin zone continues to be satisfac- 
tory. 
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(0,0) (ji/2,1/2) (it,it) <«.« n ) WO («".») «W) (0,7t) (ji/2,n/2) (u,0) 

k 

0.25 , 1 , 



(0,0) (n/2,n/2) (ji,ji) (n.»2) («,«) («n,0) (0,0) (0,u) (n/2,n/2) (it,0) 



FIG. 6. Same as Fig. 5 for t" = -J/ A. 



f s 



resulting Z ccr (k) is smalLee than the corresponding ED 
result and our SCBA oneia 

The agreement between Z ccr (k) obtained from ED and 
our SCBA approach is very good, particularly for the t— J 
model except at two k points where finite-size effects are 
evident in the ED results, and on the boundary of the 
antiferromagnetic Brillouin zone. In general, the Z ccr (k.) 
are larger if the corresponding Afc are near 
(which correspond to the ground state for all parameters 
studied here). Inclusion of other hopping terms shifts 
Afc, particularly for k ~ (jr, 0) away from \( 7r /2,w/2) > the 
corresponding weights Z c<7 (k) decrease, and the agree- 
ment between SCBA and ED results becomes poor for 
some of these wave vectors. This disagreement might 
be due to the fact that a nearest-neighbor hopping of 
the chosen sism weakens the antiferromagnetic long-range 
ordeiOEjO LJ. While the existence of this order is an es- 
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sential assumption of the SCBA, and is true for one hole 
in an infinite system, it is destroyed for a_small finite 
concentration of holes ~ 0.015 in cupratescl This con- 
centration is smaller than the corresponding one (1/32) 
in the ED studies. Another possibility for the discrep- 
ancies is that vertex corrections which seem to-not affect 
essentially the spinless hole Green function G^EJ, become 
important for Z c<T (k), particularly when the three-site 
term t" is present. 

For a quantitative comparison with photoemission ex- 
periments, it is necessary to transform the relevant op- 
erators of the appropriate multiband model, involving 
oxygen and copper orbitals, into those of the effective 
generalized t — J model. This task ffUgfyt Lp performed 
generalizing previous related studies! 41 " 



VI. SUMMARY AND DISCUSSION 
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